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We investigate the freezing phenomenon of quantum discord occurring in phase damping noise 
processes. By relating the expression of the time variation of the discord to the convex function of 
relative entropy, we obtain the necessary and sufficient conditions of the phenomenon for standard 
Bell-diagonal states. These conditions are applicable also to the phenomenon occurring in a non- 
Markovian dephasing process. Moreover, we show that the same condition and phenomenon coincide 
in a new sort of Bell-diagonal states beyond the standard form. 
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Quantum discord is regarded as a significant charac- 
teristic of quantumness [l|, y] and has attracted much 
attention in recent literatures [3|. Different from en- 
tanglement, quantum discord could capture nonclassical 
correlations arising in separable states, which has been 
found to exist in a quantum algorithm without entangle- 
ment Ij, Q . Besides the fundamental theoretic interest 
upon quantum mechanics itself, an appeal of the study 
on this subject thereby lies in the possibility to process 
quantum information in a noisy environment since pure 
entangled states are rather fragile in the real world. Re- 
cent studies on the dynamics of quantum discord in vari- 
ous noisy environment have manifested many fascinating 
features of it |6l-ll3|. It was shown that discord is more 
robust than entanglement for both Markovian and non- 
Markovian dissipative processes ^,[2]. Moreover, the dis- 
cord of a sort of Bell-diagonal states undergoing a phase 
damping process was shown to exhibit a freezing phe- 
nomenon [Ill4l3l |. i.e., the amount of discord would be 
unaffected for a finite period with the system suffering 
subsequently a sudden change from classical to quantum 
decoherence. 

A bipartite state pab is said to have nonzero quan- 
tum discord if pab can not be written as a classical- 
quantum form, i.e., pab t^ '^iPiVA){iA\ ® pf- Based 
on this benchmark, several measures have been pro- 
posed to quantify the amount of quantum discord. An 
information-theoretic measure of it [l|, y| is defined as 
Qab = I{pab) - Ic{pab), where I{pab) = S{pa) + 
S{pb) — S{pab) represents the total correlations of pab 
and S{p) = — Tr(plog2 p) specifies the von Neumann en- 
tropy; Ic{pab) = niaxn^ /[n^(pAs)] represents the clas- 
sical correlations of pab and the maximum in it is taken 
over all projective measurements 11^ on the subsystem A. 
Other quantifiers of quantum discord include the geomet- 
ric measure based on the square norm of trace distance 
[l4| , and distance- based measures either in terms of the 
relative entropy [l^ or in terms of the fidelity [Ig] . 

Generally, the dynamics of the above measures of dis- 
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cord under decoherence will manifest distinct behav- 
iors from each other quantitatively or even qualitatively, 
which suggests the complexity of the issue of decoher- 
ence from different perspectives. For the species of Bell- 
diagonal states, it is applaudable that the amount of 
discord defined via the information-theoretic measure is 
identical to that defined via the distance measure in 
terms of the relative entropy [l^- At this stage, the 
freezing and sudden-transition phenomena of discord ex- 
hibited in the Bell-diagonal states are of particular in- 
terest as the dynamical behaviors of both quantum and 
classical correlations there are consistent from different 
views related to the two kind of measures. We mention 
that comparing with those sudden-transition phenomena 
without exhibiting plateaus [SH-IO'I, the freezing dynam- 
ics of discord requires more critical conditions. Up to our 
knowledge, this freezing phenomenon was demonstrated 
mostly in systems with initial states that are of stan- 
dard Bell-diagonal form with rank two [TT| - [l3l [TtI [18| . 
A complete description of the condition related to the 
phenomenon is still lacking. 

In this paper we investigate necessary and sufficient 
conditions of the freezing phenomenon of quantum dis- 
cord occurring in the phase damping process. By elabo- 
rating the time variation of the discord of standard Bell- 
diagonal states, we associate the expression of it with 
the function of relative entropy. This enables us to give 
a complete description of conditions for the phenomenon 
by virtue of the convexity property of the relative en- 
tropy. The overall geometric structure is then depicted 
in the parametric space for the collection of states that 
satisfy these conditions. It is shown that the derived con- 
ditions are applicable also to the phenomenon occurring 
in the non-Markovian dephasing process. Finally, we dis- 
play that the same condition and phenomenon coincide 
in a new sort of Bell-diagonal states beyond the standard 
form. 

Time variation of quantum discord of Bell-diagonal 
states under phase damping. The phasing damping dy- 
namics we are going to investigate characterizes the type 
of quantum noise process that induces loss of quantum 
coherence but without energy exchange. For a two-level 
system, the dynamical evolution of the state under the 



phase damping channel is described as 



Pit) 



f pn{0) rf(i)pi2(0) 

{ d{t) P2li0) P22(0) 



(1) 



where d{t) depicts the degradation of the coherence and 
it often takes an exponential form for a Markovian de- 
phasing process. Hereafter, we shall introduce a notation 
(7=1 — d{t) and treat it a parametrized time. 

Let us consider the local phase damping process on a 
two-qubit system that initially inhabits in the standard 
Bell-diagonal state 



Pab 
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Qa,^®CTf) = ^A,|V'i)(V'i 



(2) 



where (Ti's are Pauli operators and \ipi)^s are standard 
Bell bases of the form |V'i,3> = 7f(|00) ± |11)), IV'2.4) = 
4=(|01) ± |10)). The corresponding eigenvalues A^ > 
(j = 1, • • • ,4) relate to the parameters q's via 



Ai.: 



[l±Ci=FC2+C3], A2 



-[l±Ci±C2-C3]. (3) 



As one of the subsystems Aoi B passes through the chan- 
nel of Eq. (dJ, the evolving state pAsiq) will remain to 
be of form of Eq. ^ and the corresponding eigenvalues 



Ai,3(g) = Ai,3± -[A3 - Ai], 

A2,4(g) = A2,4 ± -[A4 - A2]. 



(4) 



The amount of discord defined via the two kind of mea- 
sures, the information-theoretic measure and the distance 
measure of relative entropy, coincides for this state and 
it can be derived analytically as [1^ 

4 

QAB{q) = 1 + ^ \i{q)log2 K{q) + h2{xM)- (5) 

i=l 

Here, h2{x) accounts for a binary entropy function 
h2{x) = — xlogaX — (1 — x)log2(l — x), and XM{q) — 
[l + CM{q)]/2 in which cniq) = max{|ci(g)|, |c2(g)|, |c3|} 
and ci,2(<7) = (1 - q)ci.2- 

The variation of the discord depends on its local an- 
alytical property, recorded as QJlg{q) [i — 1,2,3) for 
the corresponding time intervals in which CM^q) equals 
to |ci(g)|, \c2iq)\, or |c3|, hence xniq) = Ai(g) + X2{q), 
A2(g) -t- Xsiq), or Ai -f A3, respectively [20]. Consider 
first the situation of CM{q) = |c3|. In this case the term 
h2{xM) in Eq. ([5]) becomes a constant, therefore one 
obtains 



^«-«" - -^''-""^ 



(6) 



It is direct to observe that for a Markovian channel, the 
above variation is always negative since the entropy of 



PAsiq) should be monotonically increasing under such a 

dephasing channel. In order to derive further the vari- 

(1 2) 
ation of Q\g (q) , we introduce the notation of relative 

entropy: Hn{x\\y) = J2'Li ^i(log2 a^i/j/i), where {xj and 

{yi} stand for two n-component random variables. By 

virtue of the following relation 
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"2(1-9) 
we are able to work out that 



[A3(<z)-Ai(<z)] 
[A4(g)-A2((7)] 



(7) 



1 



H^«-<'>=2,1-,) 



1^2(^11.-) - H,Kh)l. (8) 



where /i(i^) and ^(r/) account for probability distributions 
specified as 

p.{q) = {xM, 1 - xm}, v{q) = {1 - xa/,xm}, 

e(9) = {Ai(g),A2(g),A3(<z),A4(g)}, (9) 

,7(q) = {A3(g),A4(q),Ai(<z),A2(g)}. 

(2) 
A similar expression can be obtained for dQ\^{q)/dq 

merely through exchanging the roles of Xi{q) and A3(g) 
in Eq. §^. 

The explicit expression of the local variation of dis- 
cord derived above enables us to analyze its dynamical 
behavior conveniently. It is readily seen, from the joint 
convexity property of relative entropy [21i] , that the terms 
in Eq. ([SJ satisfy 

H2MW) < aH2{^i\\Vi) + PH2iMm) = Hmv), (10) 

where a = Xi{q) + Xsiq), 13 — X2{q) + A4(g), and £,i and 
rji are variables of probability distributions given by 



6(9) = -{Ai(q),A3(g)}, 771(g) = l{A3((7),Ai(g)}, 
a a 



6(9) = ^{A2(9),A4(9)}, 772(9) = ^{A4(9),A2(9)} 



(11) 



Thus one derives d(5^B(g)/dg < in all time intervals for 
the Markovian process. We note that this monotonicity 
of quantum discord is a natural consequence since the 
Markovian dephasing channel is always unital 116]. 

Conditions for the freezing and sudden-transtion dy- 
namics of quantum discord under phase damping. 

Theorem. The Bell-diagonal states of Eq. ^ under- 
going the phase damping process can exhibit zero quan- 
tum discord rate, i.e., the so-called freezing phenomenon 
of discord, if and only if the parameters A^ 's satisfy 

A1A4 = A2A3, (Ai - A4)(A2 - A3) > 0, (12a) 
or AiA2 = A3A4, (Ai-A2)(A4-A3) >0, (12b) 

alternatively. 



We specify that the inequaUties in Eqs. (J12ap and 
(|12bl) figure out different order relations for A^'s. For 
instance, Eq. (|12a|) indicates the possible order as: 



Ai > A2 > A3 > A4, 
A3 > A4 > A\ > A2, 



A4 > A3 > A2 > Ai, 
A2 > Aj > A4 > A3. 



(13) 



To prove the theorem, let us suppose that one starts 
from a state of Eq. ^ in which the parameters A^'s sat- 
isfy one of the order relations of Eq. (fT3|) . Note that 
any of these order relations yields that cm — |ci|, hence 

Qab = Qab^ ^'^^ ^^^ initial state. Consequently, the lo- 
cal variation of Qab (<z) during the starting time interval 
is illustrated by Eqs. (|8l)- (fTT]) . From the fact that the 
relative entropy H2{fJ.\\v) is a strict convex function, the 
equality in between the first two expressions of Eq. PH)) 
holds if and only if there are 

/^(g) = 6(9) = 6(g), iy{q) = m{q)=m{q)- (14) 

This leads clearly to the relation A1A4 = A2A3, which is 
responsible for dQAB{0)/dt = of the specified initial 
state. It is crucial to note that according to Eq. (|4]), the 
relation Xi{q)X4{q) = X2{q)X3{q) will sustain during the 
evolution. The zero-rate dynamics of discord thus will 
continue until the presumed order relation is broken. In 
addition, it is direct to observe that a similar discussion 
is applicable for the condition of Eq. (J12b[) if one starts 
from a state with ci\i — |c2|. These facts together with 
that dQ)^g{q)/dq 7^ [cf. Eq. ([6])] then complete our 
proof for the theorem. Straightforwardly, by substituting 
Eqs. (|12ap or (I12bp into Eq. ^ one obtains that 



QABiq) = 1 - h^iXi + A3) = QABiO), 



(15) 



where h2{x) stands for the aforementioned binary en- 
tropy function. 

To shed light further on the physics related to the 
above derived condition of the freezing phenomenon, we 
notice that the term H2{n{q)\\v{q)) in Eq. ^ can be 
rewritten as H2{^{q)\\i'{q)) — H4{^' {q)\\T]' (q)) , where 

,,, . ^ f XMJq) XMJq) l-XMJq) 1 -XMJq) ^ 



,, , l-XM{q) l-XM{q) XM{q) XM{q). 
V (q) = { ;; , ;; , ^r;— > ^7;— }> 



(16) 



and XM{q) = \{q) + ^2{q)- Recall that the relative en- 
tropy Hn{x\\y) tells actually how difficult it is to distin- 
guish two states of the stochastic event, {xi} and {yi}. 
The right hand of Eq. ([5]) can be understood as the re- 
duction of the distinguishability of the two states under- 
going a mixture process: £,{q) — >■ ^'{q) and rj{q) — >■ 'r]'{q), 
that is, the components of the first two of £^{q) and ri{q) 
are mixed and so well as the last two, respectively. At 
this stage, an alternative perspective on the freezing phe- 
nomenon of discord, respecting its necessary and suffi- 
cient conditions (I12ap and (|12bp . is that the distinguisha- 
bility of the two states ^(q) and 77(5) undergoing the spec- 
ified stochastic mixture process will not alter as long as 
the corresponding condition is fulfilled. 



With the necessary and sufficient conditions described 
above, we are now able to depict in Fig. 1 an overall 
geometric profile for those states in terms of parameters 
{V^i, \/A^, V^}. Regarding that < J2i=i ^» ^ 1, the 
set of standard Bell-diagonal states constitute the portion 
of a unit sphere located in the first quadrant. The surface 
responsible for the conditions in Eqs. (|12ap and (I12b[) 
comprises four leaf-shaped curves and the only common 
node of them corresponds to the maximally random state 
PAB = ^/4. The boundary of the surface figures out 
the collection of transition points which can be described 
explicitly by four independent sets of equations as 



(17) 



Ai — A2 — 


V A3 - A3, 




^3 = ^^2 = 


V Ai — Ai, 




Ai = 1 - A 


AI-A3 = (1- 


V^)' 


h = l-\ 


A^-Ai = (l- 


\^f 




FIG. 1: Geometric profile in the parametric space of standard 
Bell-diagonal states that could exhibit freezing dynamics of 
quantum discord under phase damping. 

The freezing and sudden transition phenomena of the 
dynamics are displayed schematically in Fig. 2. The tra- 
jectory of an initial state at the point / proceeds along 
the surface with a constant discord till it encounters the 
boundary at the point T, after which the discord de- 
creases to zero and the trajectory reaches finally the point 
F. At the transition point q = qx, the rate of quantum 
discord suffers a sudden change from zero to a value of 



dq I 



1 



=g+ 



2(1 -gr) 



H4{aqT)UqT)), (18) 



where we have used the relation of Eq. ([6]) and ^ and rj 
are specified in Eq. ([9]). 

It is worthwhile to point out that the conditions de- 
scribed in Eqs. (J12al) and (J12b|) are applicable also to 
the frozen discord occurring in the non-Markovian de- 
phasing process |12l |. The dynamical evolution of this 



process can still be described phenomenologically by Eq. 
(HI) , hence the dependency of the discord on the param- 
eter q does not alter. As a result, all discussions related 
to Eqs. (|12ap and (I12b[) are valid for its dynamics. Note 
that in the non-Markovian process d{t) (hence q) is not 
a monotonic function of t any more. Say, for the channel 
with random telegraph signal noise it is described as [22] 



g = 1 - d{t) = 1 - e~''*[cos(a;i) + - sin(wt)]. 



(19) 



where to = ^jAa? — ^, and a characterizes the coupling 
strength of the system with environment. Differing from 
the model concerned in Ref. Il2l in our system only one 
of the qubits is subjected to the noise channel. To dis- 
play multiple transitions and freezing phenomena of the 
discord, the parameters of the initial state herein should 
satisfy |ci(2)/c3| > e"^^^ ■ In the parametric space of 
{-yXi", VA21 V^}i we depict the evolution of the non- 
Markovian process with C1/C3 — || and 7/0; = ^jj- ^^ 
is shown that the trajectory manifests a traversing phe- 
nomenon between two separate leafs (see Fig. 2). 




FIG. 2: Illustration of the trajectory for the freezing and 
sudden transition dynamics of quantum discord in the para- 
metric space. The three blue-dashed lines account for the 
collection of Bell-diagonal states with zero discord and they 
share a common node with the two leaf-shaped curves de- 
termined by Eq. (|12a|l . The initial state is specified as 
Ai = 3/4, A2 = 3/16, A3 = 1/20. The trajectory of the solid 
line proceeding along J — >■ T — > _F describes the evolution of 
the Markovian phase damping process. For a non-Markovian 
dephasing process specified by Eq. (|19l) . the trajectory of the 
evolution proceeds further along the dotted line F — >■ T' and 
oscillates damply between the two specified leafs, displaying 
feedback and multiple transitions of the discord dynamics. 



As a final proposal of this paper, we demonstrate that 
the above described freezing dynamics of discord and as- 
sociated conditions will occur and coincide in a sort of 
new Bell-diagonal states beyond the standard form of 
Eq. (O. Consider that 



PAB 



(/4 + y^Q.af ® af) = Yl Kmm, m 



i=l 1=1 

where C12C21 7^ but all other off-diagonal elements of 
{cij} vanish. The eigenstates are shown as IV'i 3) = 
-^+|11)), |V^,4) = ;)5(|01)±e-^^-|10))with 

C12±C21 



^(100) ±< 



')± — — arctan 



C22TC11 



and the corresponding eigenval- 



ues A^'s relate to c^-'s via 



A'l 3 == T [1 ± V(cil - C22)^ + (Cl2 + C2l)2 -|- C33] 



^^2,4 ^ jt^ =*= V(cil + C22)^ + (Cl2 - C21) 



(21) 



C33J 



Suppose that the subsystem A is subjected to the local 
phase damping channel described by Eq. ^ . It turns out 
that the state pabW) will keep the form of Eq. (PH)) with 
{cy} evolving as: 033(9) = C33, and c,j{q) = (1 - g)cy 
for all other elements. Consequently, the eigenvalues of 
PAsiq), Kil)^ ^^^ derived promptly and they are verified 
to possess the similar expressions as those in Eq. (JH). 
This implies that the discord of pAsiq) should display 
the same dynamics as states of Eq. 1^. Therefore, it 
is safe to conclude that the freezing phenomenon could 
occur for the present state and the conditions of it are 
coincident with those in Eqs. (|12al) and (J12b[) . 

In summary, we have studied freezing phenomena of 
quantum discord for Bell-diagonal states under phase 
damping noise processes. We have obtained necessary 
and sufficient conditions for these intriguing phenomena 
and the overall geometric profile is then illustrated in the 
parametric space for the corresponding states. More- 
over, we have shown that the same freezing dynamics 
and conditions could occur and coincide in a sort of new 
Bell-diagonal states beyond the standard form of it. For 
further applications, our explicit derivation of the condi- 
tion for frozen discord is of potential usefulness to realize 
quantum information processing, e.g., designing quan- 
tum systems to preserve at best the quantum coherence 
in the presence of the noisy environment. 

This work was supported by the NSFC under grant 
No. 10874254. 
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